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ABSTRACT 

With the advent of large dedicated Type la supernova (SN la) surveys, K- 
corrections of SNe la and their uncertainties have become especially important in 
the determination of cosmological parameters. While A-corrections are largely 
driven by SN la broad-band colors, it is shown here that the diversity in spectral 
features of SNe la can also be important. For an individual observation, the 
statistical errors from the inhomogeneity in spectral features range from 0.01 
(where the observed and rest-frame filters are aligned) to 0.04 (where the observed 
and rest-frame filters are misaligned). To minimize the systematic errors caused 
by an assumed SN la spectral energy distribution (SED), we outline a prescription 
for deriving a mean spectral template time series which incorporates a large 
and heterogeneous sample of observed spectra. We then remove the effects of 
broad-band colors and measure the remaining uncertainties in the A-corrections 
associated with the diversity in spectral features. Finally, we present a template 
spectroscopic sequence near maximum light for further improvement on the K- 
correction estimate. A library of ~ 600 observed spectra of ~ 100 SNe la from 
heterogeneous sources is used for the analysis. 
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Introduction 



Pioneering SNIa surveys led to the surprising discovery that the cosmic expansion is 



Garnavich et al.l Il998l ; iRiess et al.l Il998l ; ISchmidt et al.lll998l ). Although much theoretical 



Perlmutter et al. 


1997, 


1998. 


1999 



work has been done on dark energy, its true nature remains a mystery ([Peebles fc Ratra 
2003T I. 



The measuremen t of the equation-of-sta te parameter can exclude certain theoretical 
mode ls of dark energy flHuterer fc Turnerll200lh . Large surveys such as the SNLs3 ( Ustier et al 



2006h and ESSENCE^ fcrisciunas et alil2005h aim to constrain the equation-of-state param- 



eter of dark energy with photometric and spectroscopic observations of hundreds of high- 
redshift SNe la. With such a large sample size, the systematic uncertainties in i^-corrections 
and photometric calibration become comparable to the intrinsic dispersion in the error bud- 
get. Improvements in the i^-correction calculations and more sophisticated estimates of the 
associated errors are thus especially important. 

A i^-correction converts an observed magnitude to that which would be observed in 
the rest frame in another filter, allowing the compariso n of the brightness of SNe la at 
various redshifts (jOke fc Sandagd Il96a ; iHogg et al.l 120021 ) . The ^-correction calculations 
require the SED of the SN la. A spectral template time series of SN la is usually used as 
an assumed SED. This is justified as there exists rema rkable homogeneity in the observed 
optical spectra of "normal" SNe la (jBranch et al.lll993l ). The dependence of i^-corrections 
on the spectral template is large at redshifts where the observed and the rest-frame filters are 
misaligned. While the observed spectra of "normal" SNe la are apparently uniform, subtle 



differences in feature streng ths and velocities do exist (INugent et al.l 1 19951 ; iBenetti et al. 



20051 ; lHachinger et al.l 120061 ) . How well the spectral template represents the SED of the 
SNe la at the misaligned redshifts therefore becomes especially important. 



Previous papers on .fT-corrections of SNe la developed met hods which made it 



possible to 



use SNe la for the determination of cosmological parameters. lHamuy et al.l (119931 ) explored 
the possibility of using SNe la as distance indicators and presented single-filter i^-corrections 
for supernovae out to redshift z = 0.5 using spectroscopic observations of SN 1990N, SN 
1991T and SN 1992A. At h i gh red shifts, considerable extrapolation is required for single-filter 
^-corrections. iKim et al.l (119961 ) developed the cross-filter K-correction to take advantage 



1 SuperNova Legacy Survey; see http://www.cfht.hawaii.edu/SNLS/ 

2 Equation of State: SupErNovae trace Cosmic Expansion; see http://www.ctio.noao.edu/essence/ 
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of the overlap at hig h redshifts of a rest-frame filter band and a redder observed filter band. 



Nugent et al.l (120021 . hereafter N02) presented a SNIa spectral template time series and a 
recipe for determining i^-corrections. The importance of broad-band colors was particularly 
emphasized. 

While if-corrections are largely driven by SN la colors, it is shown here that the diversity 
in spectral features of SNela can also be important. This paper investigates the effects of 
SNIa inhomogeneity on the determination of i^-corrections, paying close attention to the 
effects of inhomogeneity in spectral features. There are two goals for this paper. The first 
goal is to demonstrate a method for combining a large and heterogeneous set of observed 
spectra to make a representative mean spectral template. The second goal is to measure the 
remaining uncertainties in the i^-corrections once the effects of colors have been removed. 
We test the efficacy of the methods using a large library of ~ 600 observed spectra of ~ 100 
SNela. 

In Section [2], we outline the method for calculating the i^-correction of a SN la, taking 
into account the colors of the SNIa. The library of observed SNIa spectra is described in 
Section [3J In Section HJ we outline the method for color-correcting a supernova spectrum, 
a procedure necessary for fT-correction calculations and for quantifying the effects of SNIa 
spectral diversity independent of broad-band colors. Section [5] provides the prescription for 
building a mean spectral template time series which has statistically representative spectral 
features. Section [H] characterizes the associated ^-correction errors; Section [7| quantifies 
the effects of these errors on cosmology. In Section [HI we present a template spectroscopic 
sequence and its potential for further improvement on the i^-correction estimate. 



2. /^-corrections 

.fT-corrections of SNela rely on the SED of the supernova. High signal-to-noise spec- 
trophotometry of high-redshift SNela at multiple epochs is not currently feasible, and it is 
therefore necessary to use a spectral template time series. 

For a high-redshift supernova, the light of a rest-frame filter band is redshifted to a 
longer wavelength, and is observed through an observed filter band that overlaps, but is not 
i dentical to, the redshifted rest-frame filter band (Figured]). The cross-filter .fT-correction 



(IKim et al.lll996l ). K xy , allows one to transform the magnitude in the observed filter, y, to 



the magnitude in the rest-frame filter, x: 
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J\ xy (l,Z,p)- /.O lOg ^ j xTy(\)Z(X)d\ J /.X 
^Z.dlOg ^ j AT s (A)[S(A/(l+*),t^)/(l+*)]dA ) ■ 

The i^-correction, is a function of the epoch t, the redshift z, and a vector of 

parameters, p, which affect the broad-band colors of the SNIa (e.g., light-curve shape and 
reddening). Here, S designates the SED of the supernova, and A designates the wavelength. 
T x and T y denote the effective transmission of the x and y filter bands, respectively. Z denotes 
the SED for which the x — y color is precisely known. Equation [1] essentially compares the 
supernova fluxes in a rest-frame filter x and a de-redshifted observed filter y. Note that the 
normalization of the SED is arbitrary. 

In this paper, we concentrate on the .fT-corrections to rest-frame B filter band. We 
use the effectiv e transmission of four MegaCam filters, gM r Miu z Mi as examples of observed 



filters and the iBesselll (119901 ) realization of the Johnson-Morgan B filter. MegaCam is a 



wide-field imager on the Canada- France-Hawaii Telescope on Mau na Kea. The Mega Cam 



filters are closest to the US Naval Observatory (USNO) filters of ISmith et al.l (120021 ) and 
similar to the Sloan Digital Sky Survey (SDSS) filters. The methods of building a spectral 
template and estimating errors are general and not restricted to this set of filter bands. 

.fT-correction calculations also depend on the supernova redshift through the alignment 
of the rest-frame and the observed filter bands. The observed filter band progressively 
shifts toward bluer parts of the spectrum as a supernova is redshifted to progressively longer 
wavelengths (from top to bottom panel of Figure [T]) . When the observed and rest-frame filter 
bands are misaligned, as in the top and the bottom panels of Figure [TJ the i^-corrections 
are more dependent on the assumed spectral template. 

.ff -correct ions largely depend on the broad-band colors of SNela and are thus sensitive 
to anything that affects the continuum of the SED, such as Milky Way and host galaxy 
extinction, and the intrinsic colors of the supernova. When a spectral template is used as 
an assumed SED for the supernova, its continuum must be adjusted to have the same colors 
as the supernova before the X-corrections can be determined. N02 demonstrated that two 
SNe la with very different spectral features can yield similar .fT-corrections when the spectra 
of the SNe la are adjusted to the same broad-band colors. Adjusting the colors of the spectral 
template to match those of the particular SNIa significantly improves the accuracy of the 
i^-correction (N02). 

I t is an observationa l fact that SNela show a range of feature velocities and strengths 
(e.g.. iBenetti et al.ll2005l ); therefore, it will not generally be true that two SNela with iden- 



tical broad-band colors have identical SEDs. The diversities in individual feature strengths 
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and velocities between SNela must affect i^-correction calculations at some level. With 
large, modern surveys, the precision and accuracy requirements for the i^-corrections are 
only increasing. The uncertainties caused by the diversities in spectral features become 
significant in such cases. This paper will focus on the effects of spectral features, rather 
than colors, on K-corrections and address the issue by building a statistically representative 
spectral template. 



3. Library of spectra 

To study the effects of using an assumed time series of spectral templates on the deter- 
mination of .fT-corrections, it is important to use a statistically significant sample of observed 
SN la spectra which spans a wide range of supernova properties. We have constructed, and 
are continuing to update, a library of observed SNIa spectra which includes ~ 600 spectra 
from ~ 100 SNe la for this study. This paper is primarily a methods paper. Our template 
spectra, and references to the component spectra that enter into calculating the templates, 
can be found on the wetJE The spectral template time series is being constantly improved 
as more and better SN la spectra become available. The properties of the library as used in 
this paper are illustrated in Figure [2j 

The observed spectra were gathered from heterogeneous sources with a variety of wave- 
length coverage, redshifts and quality. A large fraction of the spectra was obtained from 
public resources, such as the SUSPECT0 database and publications. Other spectra were 
supp lied by past and on going supernova surveys, such as the C alan/Tololo Supernova Sur 
vev Jrlamuv et aD 1996 ). the Supernova Cosmolog y Project (Garavini et all 200 



2005; 



Lidman et aU2005r). the C arnegie Supernov a Project (IHamuy et ali2006r) . the SNLS ( 



Bronder et al.l 120071 ) and ESSENCE flMatheson et al 



2005 



Howell et al 



2005 . 



Each observed spectrum is inspected by eye for quality. Spectra which have their spectral 
features dominated by noise are discarded. The spectra taken from ground-based telescopes 
have the telluric features removed. The edges of each spectrum are inspected and trimmed 
where the flux calibrations appear unreliable. 

The library spectra are overwhelmingly of low redshift SNela (top left panel of Fig- 
ure H]). Spectroscopic observations of low-redshift supernovae are less affected by host galaxy 
contamination and are generally of higher signal-to-noise than high-redshift ones. On the 



3 See http://www.astro.uvic.ca/~hsiao/kcorr/ 
4 See http://bruford.nhn.ou.edu/~suspect/index.html 
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other hand, high-redshift surveys produce spectra with better observations of the difficult 
UV region, which is crucial for ^-correct ions to B at high redshifts. We include spectra 
from supernovae with redshifts up to 0.8, and work under the assumption that there is no 
significant evolution of spectral f eatures between low-redshift and high-re dshift supernovae 
jFolatellilbooi : lliook et alihood : iBlondin et aDhood : bronder et aDl2007h . The histograms 
in Figure [2] are separated by redshift at z = 0.2 to demonstrate the difference in character- 
istics between the two samples. 

The library spectra best cover rest-frame B and V filter bands (bottom right panel of 
Figure [2]). We will focus our analysis on i^-corrections to the B band, since it is the most 
important and commonly used band in current cosmological analysis. When the observed 
and the rest-frame B filter bands are misaligned, the spectral features in the neighboring U 
and V filter bands are especially important. The relative shortage of UV spectra is slightly 
improved by the inclusion of high-redshift spectra. 

The rest-frame epochs of the spectra are relative to B band maximum light and are 
obtained mainly from the light curves of the supernovae. When the photometry of a SN la is 
unavailable or unreliable, an estimate of the epoch is made by fitting the spectrum to other 
spectra of known ages. The distribution of the epochs peaks at maximum light (bottom left 
panel of Figure [2]). There are no late time spectra for the high-redshift sample. 



The correlation between the peak luminosity and the light-curve shape (jPhillipslll993l ) 
is a useful parameter to characterize the diversity of the library SNela. For each SNIa 
with well sampled light curves, a stretch factor is computed. The stretch factor, s, linearly 
stretches or contracts the time axis of a template light curve around the epoch rela t ive to 



1999 



maximum B band lig ht to best fit the observed light curve (jPerlmutter et al.l 11997 
Goldhaber et al.ll200ll ). The stretch factors of the library SNela span a wide range of values 
with most of the supernovae around stretch of unity (top right panel of Figure [2]). The 
high-redshift sample has a slightly higher median stretch (s = 1.049) than the low redshift 
sample (s = 0.967). This is probably caused by the preferential selection of brighter (higher 
stretch factor) SNe la at high redshifts. The median stretch factor for the combined sample 
is 0.983. 

Even t hough there are es tablished relations between light-curve shape and intrinsic SN la 
colors (e.g.. lRiess et al.lll996l ). the spectral feature strengths and velocities of SNela do not 
follow a simple relation at all ages and wavelengths. SNe la with the same light-curve shape 
do not necessarily exhibit the same spectral feature strengths; therefore, the diversity of 
spectral features cannot be completely characterized by one parameter alone. We can only 
strive for a truly diverse and representative library by including as many observed spectra 
from as many SNe la as possible. 
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4. Color-correction 

To first order, i^-corrections are driven by the diversity of SN la colors, whether this is 
attributed to extinction or to the variations of SNIa SED with age and light-curve shape. 
We illustrate the relation between /^-corrections and supernova colors in Figure [3j The 
/^-corrections are calculated for each library spectrum. The spectrum colors are determined 
from synthetic photometry on the actual observed spectrum. 

We use the colors from the spectra themselves instead of photometric colors for the 
following reasons. Many library spectra are not spectrophotometric; therefore, using colors 
from photometry introduces errors which do not originate from supernova feature inhomo- 
geneity. This choice was made also for practical reasons as the photometry of the library 
supernovae is often unavailable or unreliable. 

In Figure [3J we demonstrate the relation between the i^-correction and broad-band 
colors at two redshifts for %m i— ► B. At a redshift of 0.75 where the B and %m filter bands 
are aligned (middle panel of Figure [TJ, the i^-corrections are mostly constant with respect 
to the colors of the supernovae (top panels of Figure [3J. At z = 0.9, the filter bands are 
misaligned (bottom panel of Figure [T]), and the ^-corrections show a strong dependence on 
broad-band colors with a larger scatter around the trend (bottom panels of Figure [3]). At 
these misaligned redshifts, the .fT-correction involves more extrapolation and hence is more 
reliant on the details of the SED. The colors of the SED describe the bulk of the relation. 
The remaining scatter can be attributed to the inhomogeneity in spectral features, which is 
the focus of this paper. 



How should one allow for color v ariati ons when calculating fT-corrections? N02 utilized 

v 



the reddening law of ICardelli et ajj (119891 ) to perform the color-correction on the spectral 
temp l ate. This approach was also adopted by the High-z Supernova Search Team (IRiess et al. 



19981 : iTonrv et all 12003 ). This choice is based on the observation that SNIa color relations 
roughly follow the reddening law. The correction applies a slope as a function of wavelength 
to the spectral template in order to obtain the desired B — V color. 

Large surveys, such as the SNLS, produce well sampled multi-filter light curves and in 
turn produce multi-color information for most SNela. To utilize fully the multi-color infor- 
mation, we have developed a "mangling" (color-correction) function in lieu of the reddening 
law slope correction. Th e mangling functi on is similar to the treatment of spectral tem- 



plate color adjustment in IRiess et al.l (120041 ) . It defines splines as a function of wavelength, 



with knots located at the effective wavelengths of the filters. Non-linear least-squares fit- 
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tingj ([More et al.l 11980) is used to determine the spline which smoothly scales the spectral 



template to the correct colors. 

Figure H] shows an example of color-correction which contrasts the mangling function 
and the slope correction of the reddening law. The two methods exhibit similar behavior 
between the B and V filter bands. When the V — R color information is included, the 
scale in the R band determined by the mangling function diverges from a slope correction to 
incorporate the extra color information. The mangling function anticipates the possibility 
that the reddening law does not completely characterize all the colors of a SN la; it effectively 
pieces together segments of slope corrections smoothly to satisfy the multi-color information. 

The variations in large features of SNe la, such as the Ca H&K feature in the UV, can 
significantly alter the colors of the supernovae. It is difficult to disentangle color and feature 
variations in these cases; however, the results of the analysis in this paper are unaffected 
as long as the distinctions between the effects of color and feature are consistently defined. 
Before spectra are compared with one another, they are color-corrected to the same broad- 
band colors using the mangling function described here. This procedure consistently defines 
an effective "continuum" for all spectra and enables the study of the effect of spectral fea- 
tures on fT-corrections independent of colors, whether the color differences are caused by 
extinction, intrinsic colors or observational effects. 

Returning to Figure [31 we can see that it illustrates not only the general trend of K- 
correction in terms of supernova colors, but also the scatter around the trend which is caused 
by feature variation between supernovae. The scatter becomes more significant as the rest- 
frame filter band is redshifted away from the observed filter band (from top to bottom panels 
of Figure [3]) showing the heavier reliance on the spectra and the inhomogeneity of the spectral 
features. The goal now is to construct a spectral template time series that is representative 
of the observed SNIa spectral features and to characterize the errors associated with the 
spectral inhomogeneity observed in the scatter. 



5. Constructing a spectral template 



The spectra of normal SNe la are remarkably uniform in the time evolution of their 
spectral features ([Branch et al.lll9931 ). but some differences do exist. Some of these differences 



are associated with the var iation in light-curve shapes (INugent et al.l Il995l ; iBenetti et al. 



20051 ; lHachinger et al.l 120061 ) . while the origins of others remain unknown. 



http: / /cow. physics .wisc.edu / ~craigm / idl / idl.html 
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N02 introduced a spectral template time series that is based on one well-observed su- 
pernova spectrum for each epoch and wavelength interval. The template is constructed by 
assembling well-observed spectra in a two-dimensional grid of flux as a function of epoch 
and wavelength. The temporal gaps are filled with a simple interpolation. 

Supernova spectra with the same broad-band colors do not necessarily have the same 
spectral features inside a broad filter band. The spectral features of the template used 
for i^-correction calculations of a particular set of SNela should be representative of the 
population. The template constructed by the method outlined in N02 can cause systematic 
errors if the template represents the extremes in spectral features. In this section, we outline 
a prescription for constructing a spectral template time series of SNIa with representative 
spectral features. 

Simply averaging a large number of spectra tends to wash out spectral features, when 
there exists a sizeable range of feature velocities. Instead, we adjust the features of a base 
template to the weighted mean of the measured feature strengths from a large sample of 
observed spectra. The procedures of building the spectral template time series have been 
built into a fully automated IDL program. The program can incorporate new additions to 
the library, adjust the grid size accordingly and create a new spectral template. We outline 
the procedures in the following summary, while the details of each step are described in the 
following subsections: 

1. Assign a two-dimensional grid of epoch and wavelength bins. 

2. Color-correct each library spectrum to the same broad-band colors in each epoch bin. 

3. Measure the feature strength in each wavelength bin of each library spectrum. 

4. Assign weights to each feature strength of each library spectrum. 

5. Determine an effective epoch for each epoch bin. 

6. Determine a weighted mean feature strength for each grid element. 

7. Build a base template from a subset of library spectra. 

8. Adjust the base template spectrum to the weighted mean feature strengths using the 
mangling function. 
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5.1. Epoch bins 

We first divide the library spectra into epoch bins. When assigning an epoch bin, two 
competing factors affecting the bin sizes must be considered. First, a statistically significant 
sample size of spectra in each epoch bin is required to obtain representative characterizations 
of the spectral features; larger bin sizes yield more spectra per bin. Second, the temporal 
evolution of the spectral features within each bin should be kept as small as possible; small 
epoch bin sizes are preferred in this case. The most direct way to improve this situation 
is to increase the number of library spectra such that the epoch bin sizes can be narrowed 
without sacrificing the generality of the sample in each bin. For the current sample size, we 
cannot afford to set a bin size of one day for all epochs. Instead, we adopt variable bin sizes. 

The temporal evolution of spectral features is plotted in Figure [51 The SEDs at different 
ages are normalized to the same B band flux to emphasize the evolution of the spectral 
features. The spectrum of a SNIa evolves rapidly from the time of explosion to around 
20 days past maximum light. These epochs are also the most important when fitting light 
curves. The epoch bin sizes for these epochs are kept small such that the effects of temporal 
evolution on spectral features are minimized. The temporal evolution of feature shapes in 
the B and V bands slows down considerably beyond the age of 30 days. At these epochs, 
we adopt larger bin sizes to compensate for the small number of spectra without losing too 
much information about temporal evolution. Table [T] lists the epoch bins adopted for this 
particular set of library spectra. 

It is unclear whether the SED of a SN la with stretch s at time t relative to maximum 
light corresponds to the template of stretch s = 1 at age t or s x t. This uncertainty however 
has little effect on our template building, as the large epoch bins at late times render the 
two options indistinguishable. 

5.2. Measurements of feature strengths 

We define a set of artificial non-overlapping narrow-band pass filters for measuring 
the spectral feature strengths (Figure [H]). The bandwidths of the narrow-band filters are 
logarithmic in wavelength and are designed to be on the order of the sizes of supernova 
features. We adopt a ratio of AA/A = 0.04, which yields a bandwidth of 160A at A = 4000A 
and a bandwidth of 280A at A = 7000A. The large expansion velocities of SNe la mean that 
the bandwidth of the narrow-band filters can be set large enough such that modest noise in 
observed spectra has minimal effect on the measurements of feature strengths. 

Before the feature strengths are measured, each library spectrum is color-corrected us- 
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ing the mangling function to have the same broad-ban d colors as a supe rnova with a set 



stretch value (conventionally s — 1). The templates of iKnop et al.l (120031 ) are used to de- 
termine the colors. It is worth noting here that the procedure for calculating i^-corrections 
presented in this paper is independent of the stretch-color relation. Even though the relation 
at a particular stretch value is used here as the standard colors for the library spectra, the 
spectral template would later be color-corrected to match the colors of the particular SN la 
in question. The color-correction procedure provides a consistent "continuum" for all the 
library spectra and removes the dependence on color varying factors, such as stretch, extinc- 
tion and flux calibration. The spectral diversity of the library SNe la can now be adequately 
quantified independent of broad-band colors. 

The relative flux between neighboring features of a library spectrum is measured as 
the narrow-band colors from synthetic photometry using neighboring narrow-band filters. 
Measuring relative feature strengths eliminates the need for the absolute flux of the spectra. 
The narrow-band colors are organized in a two-dimensional grid of wavelength and epoch. 
The sizes of the grid elements are defined by narrow-band filter bandwidths and epoch bins. 
Figure [7] illustrates a schematic of the grid and specifies the number of narrow-band color 
measurements for each grid element. A typical grid element in the B and V filter bands has 
~ 20 measurements. 



5.3. Weighting feature strength measurements 



The heterogeneous nature of our library means that we need to consider the differences 
between the spectra carefully. Before the narrow-band colors are averaged in each grid 
element, some simple weighting schemes are applied to each spectrum to ensure that the 
resulting template is not dominated by peculiar spectra, supernovae with many observations 
or supernovae with extreme stretch values. Weights as a function of wavelength are also 
assigned to deal with the differences in the spectral coverage and the location of the telluric 
lines. 

We chose to include spectroscopically peculiar SNe la for the following reasons. Including 
peculiar spectra in the template and the analysis gives a more complete description of the 
population of the SNe la. Moreover, people tend to make more detailed observations of more 
peculiar SNe la. Including these spectra helps to increase the number of spectra in wavelength 
and epoch intervals which are rarely covered. To adequately include these peculiar spectra 
in the spectral template, we adopt the weight, Wt ype - For the peculiar spectra, we chose a 
relative weight of one fifth that rea sonably characterizes the int rinsic fraction of peculiar 
supernovae in the SNIa population (jBranchl l200ll ; iLi et al.ll200ll ). For a grid element with 
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Nnormai normal spectra and N pecu i iar peculiar spectra, the weight for the type of spectra is: 



For a SNIa with multiple narrow-band color measurements in a grid element, the mea- 
surements are averaged and weighted as one measurement (w mv iupiicUy) ■ 

Narrow-band color measurements at the edges of a library spectrum are assigned less 
weight than the measurements at the center. This step goes beyond the edge trimming, 
described in Section [3j in minimizing the effects of miscalibration at the edges of observed 
spectra. It also gives higher weights to library spectra which have wider wavelength coverage 
and thus have their continua better adjusted by the mangling function. The weight for the 
filter coverage, w covera ge, is a triangular or trapezoidal function in wavelength space with a 
value of unity at the effective wavelengths of the available filter coverage bands and one fifth 
at the ends of the spectrum. 

Even though the telluric features for each spectrum taken at ground-base telescopes 
are removed, we apply lower weights at the location of the telluric features to account for 
the errors associated with this procedure. The diversity of redshifts in our library means 
that most of the grid elements would have at least a few feature strength measurements 
uncontaminated by telluric features. The weight assigned for telluric line contamination, 
Wteiiuric, is inversely proportional to the equivalent width of a telluric feature model inside 
the two adjacent narrow-band filters in question. 

To ensure that a grid element is not dominated by supernovae of extreme stretch values, 
the narrow-band color measurements are weighted such that the weighted mean of the stretch 
values in each grid element is close to a representative stretch value s e ff ec ti ve (conventionally 
unity). We use a normal distribution centered on the effective stretch, s e ff ec ti ve , with a 
sigma, a s = 0.2: 



For SNela without stretch factors, the mean weight for stretch is adopted for each 
grid element. Note that we can build a template which has spectral features which are 
representative for any effective stretch value by choosing s e ff ec tive- Using spectral template 
with different effective stretch can be effective if spectral features truly vary with stretch. 
Previous studies have shown that this is the case only for certain epochs and wavelength 
regions. We will explore the topic of template spectroscopic sequence further in Section [HJ 




1 / lN norma i for normal SNela 
l/5N pecu i iar for peculiar SNela 



(SSN — Seffective) 

2^ 



21 



W stretch = exp 



13 



An effective epoch, t e f ■fective is then determined for each epoch bin using the epochs of 
library spectra, £j, and their weights, w\: 



t 



effective 



EM 



The letter % denotes the ith spectrum in the epoch bin in question. The weight, w\ 
is the product of the weights for supernova type, multiplicity, wavelength coverage, telluric 
features and stretch. The resulting effective epochs for our current library are listed in 
Table [H Within an epoch bin, library spectra with epochs tsN progressively closer to the 
effective epoch are then assigned progressively higher weights. We use a normal distribution 
centered on the effective epoch, t e ff ective , with a sigma, a t , corresponding to the size of the 
epoch bin: 



w, 



epoch exp 



(tsN — teffective) 

2^ 



The total weight for the ith library spectrum can then written as: 



^i ^type,i X W 'multiplicity ,i X Wcoverage,i ^telluric,i X Wstretch,i X W e p 0c fo i 

Assuming that the spectral features of SNe la depend on the stretch parameter at some 
level, weighting feature strength measurements in terms of epoch and stretch helps us to 
disentangle the two effects on spectral features in an epoch bin of finite size. 



5.4. Building the spectral template 

The weighted mean of the narrow-band color measurements, c, is determined for each 
grid element (in epoch and in wavelength) using the total weight, Wi, and the narrow-band 
color measurements Cj for the i spectrum: 

c _ Ej w i°i 

EiWi 



The grid of weighted mean narrow-band colors in epoch and wavelength now character- 
izes the spectral features of the desired SN la spectral template time series. 
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It is important to start with a base template spectrum with the correct feature shapes. 
We build a base template spectrum at the effective epochs using a selected set of high 
signal-to-noise spectra from the library. The coverage in epoch and wavelength is kept below 
three spectra per day per wavelength to avoid washing out the spectral features. Each library 
spectrum is corrected to the same colors and the same B magnitude at each epoch. Gaussian 
weights centered on the effective epoch are applied to the flux. The weighted mean flux is 
then computed for each effective epoch. 

The grid of weighted mean narrow-band colors, c, as a function of wavelength and epoch 
is now used as the input color information for the mangling function described in Section HI 
Each spectral feature of the base spectrum is adjusted by the mangling function to have 
the weighted mean strength of the library spectra. The adjusted template at the effective 
epochs is then smoothly interpolated to an integer epoch grid from t = —19 to t = 85 at 
each wavelength. Finally, the spectral template is color adjusted to a light-curve template. 

Note that adjusting the flux and the colors of the spectral template to a particular set 
of template light curves is optional for the following reasons. The constant of normalization 
for the template SED is cancelled in Equation [TJ Furthermore, the broad-band colors of the 
spectral template is later adjusted to match the colors of the SNIa in question. Here, we 
choose to adjust our spectral template to a light-curve template to allow for a more natural 
interpolation of flux between epochs. 



In this section, we characterize the errors associated with the use of an assumed SNIa 
spectral template time series when calculating the i^-corrections of SNela. /^-correction 
errors are quantified here as the differences between the i^-correction calculated using an 
actual supernova spectrum S s and that calculated using an assumed template spectrum 
S T , where the continuum of S T is adjusted with the mangling function to have the same 
broad-band colors as S s (Section Hj). From Equation [H we obtain: 



This definition of i^-correction errors removes the dependence on the zero points and 
focuses on the effects of an assumed spectral template and the alignment of the rest-frame 
and observed filters. The first term in Equation [2] provides the normalization such that the 



6. Quantifying /^-correction errors 




(2) 



- 15 - 



two spectra have the same flux in the rest-frame filter band. The second term then yields 
the magnitude difference caused by the spectral feature differences between the supernova 
spectrum and the template. 

Here, we utilize the library spectra again to quantify the i^-correction errors. The library 
spectra are artificially redshifted from z — to z — 1.2 to determine the .fT-correction as 
a function of redshift. The ^-correction differences between the template and each library 
spectrum are then determined for all redshifts with various combinations of rest-frame and 
observed filter bands. We then define the statistical error as the .fT-correction differences 
added in quadrature at each redshift and the systematic error as the average of the K- 
correction differences at each redshift. 



When a relation between light-curve shape and color (e.g.. iRiess et al.lll996l ) is used, 
one has available the predictions of multiple rest-frame color information from the fit to 
the observed light curve. The spectral template can be color-corrected using the predicted 
rest-frame colors in this case. On the other hand, when multi-color observations of a SNIa 
are available, one may choose to use the more empirical approach of correcting the spectral 
template directly to the observed colors. We explore the .fT-correction errors from both 
color-correction options in the following subsections. 



6.1. Correcting the template using model rest-frame colors 

We first look at the case where the colors of the spectral template are corrected by the 
rest-frame colors. When determining the .fT-correction to a rest-frame filter, the neighbor- 
ing spectral regions become important when the rest-frame and observed filter bands are 
misaligned. To determine the errors caused by the inhomogeneity in spectral features, it is 
important to use spectra which have adequate wavelength coverage and can be accurately 
color-corrected by the mangling function. We took all the library spectra which have U, B 
and V filter band coverage (51 spectra, about 8% of our current library) to characterize the 
errors in i^-correction for the rest-frame B filter band. The selected library spectra range in 
epoch from t = — 15 to t = 65. The rest-frame colors of the library spectra are measured by 
synthetic photometry. Before the .fT-correction difference between a library spectrum and 
the template is calculated, the template spectrum is color-corrected to have the same U — B 
and B — V colors as the library spectrum. 

We randomly select half of the library spectra with U, B and V band coverage to calcu- 
late the .fC-corrections. We use the other half and the rest of the spectra without complete 
U, B and V coverage to make the spectral template, such that the spectral template is 
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always independent of the library spectra used to examine the i^-correction differences. The 
process of random selection, building a template with one half of the sample and calculating 
the iT-correction errors with the other half was repeated 100 times. 

We plot the statistical and systematic errors for .fT-corrections from observed c/mT m^m^m 
to rest-frame B in Figure [HJ The solid black curve is from the new template building method 
presented here; the dotted curve is from the revised template of N02. The thickness of the 
black curve shows the dispersion from the random selection process and the thin white curve 
is the mean. We assume here that the revised template of N02 is independent of the 51 
spectra used for the error estimate. 

Note that the errors presented here are for individual observations. Combining multiple 
data points to derive overall light-curve parameters should reduce these errors (Section [7]). 
The statistical errors are a function of redshift. They follow the pattern of a minimum at 
the redshift where the observed and the rest-frame filter bands are aligned and two maxima 
where the two filter bands are misaligned. The pattern repeats as a function of redshift 
and as one observed filter is switched to another. At a redshift where the filter bands are 
misaligned, the i^-correction error for a single observation on a light curve is important at 
the 0.04 mag level. Because the spectra are compared with consistent broad-band colors, 
we argue that the 0.04 mag error is mostly due to the inhomogeneity in spectral features of 
SNela. 

The right panel of Figure [H] shows that the new template causes lower systematic errors 
than the revised template of N02. This is because the new template has more representative 
spectral features for the 51 spectra considered here. This implies that the N02 template 
is only representative for the handful of spectra used to make it. The method presented in 
Section fallows for the inclusion of a large sample of spectra for building a spectral template. 
Assuming that our current library is a more representative sample than the sample used by 
N02, then the new template will reduce the systematic errors caused by spectral features. 

Note that we do not expect the mean systematic error to be identically zero in these 
plots. An ideal spectral template for this particular population of spectra would yield dis- 
tributions centered on zero at all redshifts and epochs. The deviation from zero for the new 
template reflects the fact that only 51 out of the ~ 600 library spectra were used (due to the 
wavelength coverage constraint) to estimate the errors and the requirement that the sample 
used to build the template is disjoint from that used to perform the test. They are consistent 
with zero within the errors. 

When a SNIa lies at a redshift where the observed and rest-frame filter bands are 
misaligned, all the observations on the light curve would have larger i^-correction errors. 
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The larger errors are attributed to the larger extrapolation and the heavier reliance on 
the spectral template. Multiple epochs may reduce the errors, but the redshift dependent 
nature remains. SNela at different redshifts would have different i^-correction errors. This 
redshift-dependent feature of the iCcorrection error must be considered in the determination 
of cosmological parameters. 

6.2. Correcting the template using observed colors 

Here we perform the same error analysis as Section 16.11 but the template spectra are 
color-corrected using "observed" colors. The observed colors are obtained from the synthetic 
photometry using deredshifted Qm^m^m observed filter bands. As the redshift of a super- 
nova spectrum increases, the observed filter bands cover progressively bluer regions of the 
spectrum. The locations of the spline knots for the mangling function are not fixed as in the 
previous subsection, but instead depend on the redshift of the supernova. 

Because the observed filter bands are shifted as a function of redshift, spectra with wide 
wavelength coverage are required for the analysis. There are 12 spectra with adequate wave- 
length coverage from seven SNela, SNe 1981B, 1989B, 1990N, 1990T, 1991T, 1991bg, 1992A, 
covering ages from —14 to 37 days relative to maximum light. The spectral template built 
for this analysis excludes these 12 spectra, such that the spectral template is independent 
of the these spectra used for the error estimate. We consider the example of .^-correction 
from observed tm band to rest-frame B band in Figure [H 

The left panel plots the comparison of statistical fCcorrection errors calculated using the 
revised template of N02 and the new template. The result is similar to that of Figure [HI Since 
the spline knots are placed at different locations for the two methods of color-correction, the 
similarity between the results gives us confidence that the method of constructing spectral 
templates presented in Section [5] is largely insensitive of the way the spectra are color- 
corrected. 

Color-correcting using observed colors has the advantage of being independent of an 
assumed stretch-color relation, but the .^-correction errors depend on the availability of 
observed colors. The effects of the availability of observed colors are shown in the right 
panel of Figure O We consider three cases where the new template is color-corrected by 
both observed gu — tm and r M — i M colors, by observed gu — ?"m color and by observed 
tm — iM color. They illustrate the point that in the treatment of iCcorrection, the errors 
should also reflect the availability of observed colors. For example, the .^-correction error 
is increased considerably for a SNIa missing g M — r M observed color at z — 0.36 (dashed 
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curve of the right panel of Figure [9]). A similar effect is also observed for a SNIa at z — 0.54 
missing r M — i M (dotted curve of the right panel of Figure [9]). 

At redshifts where the observed and rest-frame filter bands are misaligned, the observed 
color which straddles the rest-frame filter band is the most important. This is illustrated 
in Figure [101 For a SNIa that is at a misaligned redshift, and missing this observed color 
information, the associated A-correction error arises not only from inhomogeneity of SNe la 
spectral features, but also from the missing color information. 



6.3. Determining the weights from the minimization of errors 



In Section 15.31 we selected two constants for the weighting scheme, constants for w type 
and w coverage, in order to ensure that the spectral template produced is not dominated by an 
extreme of the SNIa population. We used a constant of 0.2 for w tyP e to simulate the fraction 
of peculiar supernovae in the SNIa population. We used a constant of 0.2 for w coverage to 
weight down the edges of the spectra. Alternatively, these constants can be quantitatively 
determined by minimizing the systematic A-correction errors of a set of spectra. 

We set the constants for w tyP e and w coverage as parameters which are free to vary be- 
tween and 1. We took all of the spectra which have enough wavelength coverage to allow 
for accurate mangling to determine the systematic errors in redshift as discussed in Sec- 
tion 16.11 The parameters are quantitatively det ermined by minimizing these errors using 



non- linear least-squares fitting (IMore et al.l 119801 ). At each iteration, a spectral template is 



produced from the entire spectral library with the specified parameters, and the averages of 
the A-correction differences between this template and the library spectra are taken as the 
systematic errors in redshift to be minimized. 

By minimizing the systematic errors, the optimal constants for w type and w coverage were 
determined to be 0.244 and 0.498, respectively. The constant for w type roughly reflects the 
fraction of peculiar spectra in the sample. The larger constant for w coverage reflects the fact 
that the edges of these library spectra are used to determine the errors at redshifts where 
the filters are misaligned. 

The spectral template produced by using the quantitatively determined parameters 
shows very little difference from the spectral template presented in Section [5J The templates 
at maximum light are compared in Figure [TTJ The largest A-correction (to B) difference 
between the two templates is 0.004. Even when the two constants are set as their maximum 
value of unity, the maximum A-correction difference is still less than 0.01. The spectral 
template building process is therefore largely insensitive to the way these constants are 
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determined. 

The statistical error estimates remain on the order of 0.01 at the redshifts where the 
filters are aligned and 0.04 at the redshifts where the filters are misaligned. The minimization 
of the systematic error represents an improvement of ~ 0.0005 at the misaligned redshifts 
and has no effects on the shape of the systematic errors in redshift. The insensitivity to 
these parameters is mostly due to the fact that our spectral library is large and encompasses 
SNela which span a wide range of properties. 

Despite the fact that the template building process is largely unaffected by the choice 
of these parameters, the weights determined from the minimization of errors depend on the 
characteristics of the sample used for the error estimates. To use the method presented 
in this section, one should ensure that the sample of spectra has good spectral and epoch 
coverage and is reasonably representative in properties such as the stretch and the fraction 
of peculiar spectra. Using an unrepresentative sample of spectra has the potential to skew 
the template toward a particular extreme of the SN la population. 

7. Impact of ^-correction errors on cosmology 

In the preceding sections we have concentrated on determining the iT-correction errors 
for individual observations. Since one generally observes a supernova at multiple epochs, 
it is natural to ask what the final effect of these errors will be on the derived light-curve 
parameters, such as the peak magnitude in B. 

This question is currently difficult to answer. The critical issue is how correlated the in- 
correct ion error is at different epochs. If the errors were completely uncorrelated, as is often 
assumed in the literature, then the final error would be approximately reduced by l/y/N, 
where N is the number of observations (this is not quite true because different observations 
tend to have different weights in a light-curve fit). Alternatively, if the iT-correction errors 
were perfectly correlated, then there would be no reduction in the final error from multiple 
observations. 

In reality, we are somewhere between these extremes. It is easy to understand why 
observations at different epochs would be correlated. For example, imagine a supernova at 
a redshift where a particular feature is barely outside a given observed filter for the fiducial 
SED. If this supernova has unusually high expansion velocities, then the feature might be 
shifted just inside the observed filter, which would affect the i^-correction. Since a supernova 
with a high expansion velocity relative to the fiducial template is likely to also have a high 
relative expansion velocity at a later date, the same effect will occur when the supernova is 
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observed several days later. Thus, the errors in the .^-correct ion will be correlated between 
different observations. 

What is needed to determine the level of correlation is SNela with spectra obtained 
at many epochs. One can then calculate the .fT-correction errors for the same supernova at 
different epochs and form the covariance matrix. Unfortunately, at the current time there are 
not nearly enough SNe la with multiple published spectra. N ew data samples su ch as those 



which should be provided by th e Carnegie Sup ernova Project ( Hamuy et aL 12006), Katzman 



Auto matic Imaging Telescope (jLi et al.l l2000l ). the Supernova Factory ( jWood-Vasey et al. 



2004 ) and the European Supernova Collaboration as part of a European Research Training 



Network (e.g.. iPastorello et al.l 120071 ) will help address this issue. 



Nonetheless, we can try to develop a feel for the importance of this effect with the 
current data, although this will necessarily be highly speculative. We took all of the SNe la 
with multiple spectra which have enough wavelength coverage to allow for accurate mangling, 
calculated their .fT-correction errors as discussed in Section I6.1[ and formed the covariance 
matrix by taking the appropriate products. This covariance matrix has 275 unique off- 
diagonal elements (from 25 epoch bins defined in Tabled]). The current data set provides 
information for only 10% of these, many of which have entries from only one supernova. 

Clearly, with only 10% of the elements specified, we cannot estimate the final error 
without making some additional coefficients. The distribution of correlation coefficients for 
the .^-correct ion from iu to B at z = 0.6 is shown in Figure [T21 plotted against epoch 
difference between the two bins. Note that some of the correlation coefficients are larger 
than 1, which is non-physical and is related to the fact that we only have one supernova in 
that bin from which to estimate the covariances. There are no clear trends visible, which is 
not surprising given the paucity of data. 

Therefore, to proceed, and for lack of any better model, we measure the mean correlation 
coefficient (~ 0.5) and use that to calculate all of the off-diagonal elements. We carry this 
operation out at three redshifts (z = 0.6,0.7,0.8), in all cases calculating the error in the 
i^-correction from %m to B. The alignment between B and %m is best at z — 0.8. In order 
to make this estimate, we have to have some idea how iu band observations are typically 
distributed at these redshifts. For this we use three actual SNLS SNela, 03D4dy (z = 0.604), 
04Dlsi (z = 0.701), and 04Dlks (z = 0.798). At higher redshifts, we actually tend to have 
slightly more measurements because of time dilation. 

We then calculate the expected error on the peak B magnitude in three different ways: 
assuming that the errors are completely uncorrelated, using only the covariance matrix el- 
ements we have actually measured (i.e., with 10% of them non-zero), and using the above 
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model for the unknown coefficients. The results are summarized in Table [2j Roughly speak- 
ing, it appears that the errors on the final magnitude are a factor of 2 or so larger than 
one would expect in the absence of correlations. We must emphasize that these numbers 
are tentative, but they should be a considerable improvement over the assumption that the 
K-correction errors at different epochs are uncorrelated. 

Differences between local and high-redshift SNe la are also of interest in cosmology. The 
detail ed comparison of spectra of nearby and distant SNe la is the subject of many papers 



e.g., iBlondin et al.ll2006l ; iBronder et al.l 120071 ). In the previous sections, we have assumed 



that the effect of the evolution in redshift is negligible compared to the intrinsic variations in 
spectral features between SNe la at the same redshift. Here, we perform a basic test on the 
impact of including both low-redshift (z < 0.2) and high-redshift samples in the template 
building process. We build the spectral template using only low-redshift spectra to compare 
with our new spectral template built using the full sample of library spectra (Figure [T3|) . 
The difference in .fT-corrections between the two spectral templates is small (typically 0.01 
to rest-frame B band) and mostly reflects the difference in the wavelength coverage of the 
two samples rather than true evolution. 



8. Template spectroscopic sequence 



It is possible to establish a one-to-one relation betw e en light-curve shape a nd spectral 
feature strengths (jNugent et al.lll995l ; iBenetti et al.ll2005l ; lHachinger et al.l 120061 ) for certain 
features around the epoch of maximum B band light. Quantifying the relative feature 
strengths of a large sample of SNIa spectra (Section 15.21) allows us to explore the possibility 
of constructing a template spectroscopic sequence as a function of the light-curve shape. 

As a first test, we consider the well known spectroscopic sequence of S i II A5 972 to Si II 
A6355 ratio around B band maximum light first explored by lNugent et al.l (119951 ). We select 
41 library spectra near maximum light (one day before to one day past maximum light) with 
adequate V band coverage. The spectra are of 28 SNe la with available light curves and 
stretch factors. 

We perform a principal component analysis (PCA) on the feature strength measure- 
ments to characterize the spectral features in this region. PCA reduces multidimensional 
data to principal components with fewer dimensions by identifying patterns in the data. 
PCA has been employed for analyzing astronomical spectra in a wide yariety of application s 



Ronen et al 



19991 : iKong fc Chend 1200 ll ; Madgwick et al.l 120031 : iFerreras et al.l l20Q6h . 



(e.g., 

Here, PCA is applied to our color-independent feature strength measurements (Section 15.21) 
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of the selected spectra. 

We consider the spectral region in the V band (4600 A to 6600A) . Each spectrum is char- 
acterized by nine feature strength measurements. Using the feature strength measurements 
instead of the observed spectra themselves minimizes the effects of noise in the analysis. 
The first principal component is taken along the direction of maximum variance in the nine- 
dimensional space and constitutes ~ 60% of the total variance. We assume here that the 
first principal component describes the pattern instigated by the variation in light-curve 
shapes (e.g., stretch). We establish a one-to-one relation between the pattern described by 
the first principal component and the stretch factors of the SNe la and derive a spectroscopic 
sequence. 

The resulting template spectroscopic sequence as a function of wavelength and stretch 
factor is plotted in Figure [EH The color dependence of SN la spectra on their stretch factors 
is removed by color-correcting all the spectra to the same broad-band colors. Figure [T41 then 
exhibits the spectral variations due to the varying stretch factor independent of colors. The 
feature strength ratio of Si II A5972 to Si II A6355 decreases and the equivalent width of 
the Fe trough near 4800A trough stays roughly constant as the stretch factor increases (as 
the declini ng rate decreases ). This description of the spectroscopic sequence agrees with the 



findings of lHachinger et al.l (120061 ). It represents an approximately 0.005 mag reduction in 
error for a /^-correction from %m to V at a misaligned redshift of z = 0.5 for a SN la with an 
extreme stretch value. The improvement is small compared to the 0.04 mag statistical error 
caused by feature differences which cannot be predicted by light-curve shapes. 

A spectroscopic sequence as a function of stretch factor is possible only for a limited 
range of epochs and spectral regions. SNe la with the same stretch factor and intrinsic col- 
ors in general do not have identical spectral features at all wavelengths and at all epochs. 
Some SNe la with extreme stretch values sometimes exhibit "normal" spectral features, while 
other S Ne la with moderate stretch values sometimes are identified as spectroscopically pe- 



culiar (jGaravini et al.ll2004j ). The spectral features of SNe la clearly depend on more pa- 
rameters than the stretch factor alone. This again emphasizes the need to characterize the 
i^-correction errors with as large a sample of supernova spectra as possible. 



9. Conclusions 

In the era of large SNe la magnitude-redshift surveys, the determination of fT-corrections 
and their errors becomes especially important. While i^-corrections are largely driven by 
the colors of a SN la, it is shown here that spectral features also play a significant role. The 
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statistical errors caused by the diversity in spectral features are measured to be 0.04 mag 
for a single observation at a redshift where large extrapolation is required. 

We have presented in this paper improvements in the .fC-correction calculations and in 
the error estimates. The method for color-correction simultaneously allows for the adjust- 
ment of multiple colors and a consistent continuum for the comparison of spectral features 
between spectra. The method of deriving a representative spectral template allows the inclu- 
sion of a large sample of observed spectra and reduces the systematic errors of i^-corrections. 
The redshift-dependent nature of /^-correction errors derived here shows that SNe la at differ- 
ent redshifts should be weighted differently in the determination of cosmological parameters. 

This paper is based in part on data gathered with the 6.5 meter Magellan Telescopes 
located at Las Campanas Observatory, Chile. It is also based on observations obtained at the 
Gemini Observatory and the W. M. Keck Observatory on Mauna Kea. The authors would 
like to recognize the very significant cultural role and reverence that the summit of Mauna 
Kea has within the indigenous community of Hawai'i. We are grateful for our opportunity to 
conduct observations from this mountain. Canadian authors acknowledge support from the 
Natural Sciences and Engineering Research Council of Canada and the Canadian Institute for 
Advanced Research. This research has made use of the NASA/IPAC Extragalactic Database 
which is operated by the Jet Propulsion Laboratory, California Institute of Technology, 
under contract with the National Aeronautics and Space Administration. This work is also 
based in part on Carnegie Supernova Project observations supported by the National Science 
Foundation under Grant No. 0306969. 
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Fig. I. — The pairing of observed filter % M and rest-frame filter B at redshifts 0.6 (top 
panel), 0.75 (middle panel) and 0.9 (bottom panel). The solid curves are the transmission of 
B band, and the dotted curves are the transmission of the deredshifted %m band. The filter 
bands are misaligned at z — 0.6 and z = 0.9. The .fT-corrections at these redshifts depend 
heavily on the assumed spectral template. 
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Fig. 2. — Characteristics of our current spectral library. The histograms are separated by 
redshift at z — 0.2 to demonstrate the difference in characteristics between the two samples 
(grey for z < 0.2 and black for z > 0.2). The top left panel plots the histogram of the 
redshifts of the library supernovae. The top right panel plots the histogram of the stretch 
factors of the library supernovae. Only supernovae with reliable photometry are included 
in this plot. The median stretch value for the combined sample is marked with a dotted 
vertical line. The bottom left and right panels plot the histograms of the rest-frame epoch 
relative to B band maximum light and the rest-frame wavelength coverage of the library 
spectra, respectively. 
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Fig. 3. — i^-corrections from observed %m band to rest-frame B band as a function of broad- 
band colors, U — B (left panels) and B — V (right panels). The i^-correction-color relations 
are plotted at the redshifts z = 0.75 (top panels) and z = 0.9 (bottom panels). The %u 
and B filter bands are aligned at z = 0.75 and misaligned at z = 0.9. All the library 
spectra with adequate wavelength coverage are included. Each point represents one library 
spectrum. Filled and open circles represent normal and peculiar SNela, respectively. The 
diversity in colors mostly reflects the time evolution of supernova colors. The scatter reflects 
the differences in spectral features. 
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Fig. 4. — An example of color-correction which contrasts the two color-correction meth- 
ods: reddening law slope correction and mangling function. The spectral template is color- 
corrected to match the colors of an observed spectrum. The top panel shows three spectra: 
the observed spectrum (grey), the spectral template corrected by the mangling function 
(black solid) and the spectral template corrected by the reddening law (black dotted). The 
bottom panel shows the color-correction scale for the mangling function (black solid) and 
for the reddening law (black dotted). The black dots locate the spline knots for the man- 
gling function. The mangling function adjusts the continuum using all the color information 
supplied, B — V and V — R. 
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Fig. 5. — An illustration of the time evolution of spectral features of SNela. A spectral 
template is plotted from t = — 15 to t = 48 relative to maximum B band light. Template 
spectra at different epochs are normalized to the same B band flux to emphasize the evolution 
of the spectral features. Spectral features of SNe la evolve rapidly around maximum light and 
slow down past t = 30. This emphasizes the importance of small epoch bins near maximum 
light. 



-31 - 




2000 4000 6000 8000 10000 

wavelength (A) 

Fig. 6. — The defined narrow-band filters for the measurement of spectral feature strengths. 
The filters are defined to have logarithmic bandwidths with a ratio of AA/A = 0.04. The 
bandwidths are on the order of the sizes of SN la spectral features. 
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Fig. 7. — A schematic of the two-dimensional grid defined by the epoch bins and the narrow- 
band filters. The number of narrow-band color measurements in each grid element is plotted 
in the third dimension. 
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Table 1. The assignment of epoch bins 



Lower epoch Upper epoch Bin size Effective epoch Number of spectra 
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Fig. 8. — ^-correction errors to the B band where the template spectrum is color-corrected 
with rest-frame colors. The errors of the new template and the revised template of N02 
are represented with solid black and dotted curves, respectively. The statistical errors are 
plotted in the left panel, and the systematic errors in the right panel. The thickness of the 
black solid curve shows the dispersion from the random selection process and the thin white 
curve is the mean. The errors presented here are for individual observations on a light curve. 



Table 2. i^-correction errors in peak magnitudes 



z = 0.6 z = 0.7 z = 0.8 



Uncorrelated 0.011 0.011 0.004 

Partially correlated 0.013 0.011 0.006 
Correlations model 0.023 0.018 0.010 



Note. — The approximate effects of including correlations between .fT-corrections on dif- 
ferent epochs on the peak magnitude in B. The uncorrelated numbers reflect the assumption 
that the i^-correction errors are completely uncorrelated, which is unlikely to be true as ex- 
plained in the text. The partially correlated model only includes those correlations we can 
estimate from current data (~ 10%), and the correlations model reflects the effect of very 
roughly estimating the unknown correlations using a simple model (setting them all equal 
to the mean value we can estimate). 
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Fig. 9. — Statistical errors of ^-correction from observed %m band to rest-frame B band 
where the template spectrum is color-corrected with "observed" colors. The observed colors 
are from the synthetic photometry using the observed filter bands Qm^m^m- The left panel 
compares the errors of the revised template of N02 and the new template. Both templates 
are color-corrected by observed qm — i>M and %m — tm colors before the calculation. The right 
panel shows the effect of the availability of observed colors on the errors. The solid, dotted 
and dashed curves represent the errors of the new template corrected by both Qm ~ r M and 
r M ~ iu colors, qm — f'M color alone and r M — i M color alone, respectively. 
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Fig. 10. — The location of the observed filters at redshifts z = 0.36 (top panel) and z = 0.54 
(bottom panel). We consider the same example of .fT-correction as in Figure The solid 
curves represent the rest-frame B band, while the dotted curves represent the observed 
9m^m^m bands. At redshifts where the observed and rest-frame filter bands are misaligned, 
the observed color which straddle the rest-frame filter band (highlighted for each plot) is the 
most important. 
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Fig. 11. — Comparison between the spectral template built using weights determined by 
minimizing the .fT-correction errors and preselected weights. The templates at maximum 
light are shown. The template building process is largely insensitive to the way these con- 
stants are determined. The largest .fT-correction (to B) difference between the two templates 
is 0.004. 
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Fig. 12. — Estimated correlation coefficients for the ~ 10% of the off-diagonal elements of the 
^-correction covariance matrix that we can currently measure plotted against the difference 
between the two epochs. Here the z = 0.8 %m >— ► B is shown. At this redshift the filters line 
up fairly well. Note that some of the correlation coefficients are greater than one, which is 
unphysical, and results from the fact that many of these bins have only one SNIa in them. 
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Fig. 13. — Comparison between the spectral template built with the full sample of library 
spectra (solid curve) and the spectral template built with only low-redshift spectra (dotted 
curve). The templates at maximum light are shown. The difference in i^-corrections between 
the two spectral templates is small and mostly reflects the difference in the wavelength 
coverage of the two samples rather than true evolution. 
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Fig. 14. — The template spectroscopic sequence in the V band at maximum light. The 
template is derived from 41 library spectra of 28 SNe la with a wide range of stretch factors. 
The pattern is prescribed by the principal component analysis of the narrow-band color 
measurements of the library spectra. 



